Theoretical insights into the HO●-induced oxidation of chlorpyrifos pesticide: Mechanism, kinetics, ecotoxicity, and cholinesterase inhibition of degradants.

The oxidation of the common pesticide chlorpyrifos (CPF) initiated by HO● radical and the risks of its degradation products were studied in the gaseous and aqueous phases via computational approaches. Oxidation mechanisms were investigated, including H-, Cl-, CH3- abstraction, HO●-addition, and single electron transfer. In both phases, HO●-addition at the C of the pyridyl ring is the most energetically favorable and spontaneous reaction, followed by H-abstraction reactions at methylene groups (i.e., at H19/H21 in the gas phase and H22/H28 in water). In contrast, other abstractions and electron transfer reactions are unfavorable. However, regarding the kinetics, the significant contribution to the oxidation of CPF is made from H-abstraction channels, mostly at the hydrogens of the methylene groups. CPF can be decomposed in a short time (5-8 h) in the gas phase, and it is more persistent in natural water with a lifetime between 24 days and 66 years, depending on the temperature and HO● concentration. Subsequent oxidation of the essential radical products with other oxidizing reagents, i.e., HO●, NO2●, NO●, and 3O2, gave primary neutral products P1-P15. Acute and chronic toxicity calculations estimate very toxic levels for CPF and two degradation products, P7w and P12w, in aquatic systems. The neurotoxicity of these products was investigated by docking and molecular dynamics. P7w and P12w show the most significant binding scores with acetylcholinesterases, while P8w and P13w are with butyrylcholinesterase enzyme. Finally, molecular dynamics illustrate stable interactions between CPF degradants and cholinesterase enzyme over a 100 ns time frame and determine P7w as the riskiest degradant to the neural developmental system.

New insight into environmental oxidation of phosmet insecticide initiated by HO˙ radicals in gas and water - a theoretical study.

Phosmet is an organophosphorus insecticide widely used in agriculture to control a range of insects; recently, it was banned by the European Union in 2022 due to its harmful effects. However, its environmental degradation and fate have not yet been evident. Thus, phosmet oxidation by HO˙ radicals was theoretically studied in this work using the DFT approach at the M06-2X/6-311++G(3df,3pd)//M06-2X/6-31+G(d,p) level of theory. Three different mechanisms were considered, including formal hydrogen transfer (FHT), radical adduct formation (RAF), and single electron transfer (SET). The mechanisms, kinetics, and lifetime were studied in the gas and aqueous phases, in addition to its ecotoxicity evaluation. The results show that FHT reactions were dominant in the gas phase, while RAF was more favourable in the aqueous phase at 298 K, while SET was negligible. The branching ratio indicated that H-abstractions at the methyl and the methylene groups were the most predominant, while the most favourable HO˙-addition was observed at the phosphorus atom of the dithiophosphate group. The overall rate constant values varied from 1.2 × 109 (at 283 K) to 1.40 × 109 M-1 s-1 (at 323 K) in the aqueous phase and from 6.29 × 1010 (at 253 K) to 1.32 × 1010 M-1 s-1 (at 323 K) in the gas phase. The atmospheric lifetime of phosmet is about 6 hours at 287 K, while it can persist from a few seconds to several years depending on the temperature and [HO˙] concentration in the aqueous environment. The QSAR-based ecotoxicity evaluation indicates that phosmet and its degradation products are all dangerous to aquatic organisms, although the products are less toxic than phosmet. However, they are generally developmental toxicants and mutagenicity-negative compounds.

Hydroxyl radical-initiated decomposition of metazachlor herbicide in the gaseous and aqueous phases: Mechanism, kinetics, and toxicity evaluation.

The oxidation of widely-used herbicide metazachlor (MTZ) by hydroxyl radical (HO•) in the gas and the aqueous phases was investigated in terms of mechanistic and kinetic behaviors using the M06-2X/6-311++G (3df, 3pd)//M06-2X/6-31 + G (d,p) level of theory over the temperature range 250-400 K. The formal hydrogen transfer, HO•-addition, and single electron transfer mechanisms were considered. The overall rate constants in the gas phase range from 8.40 × 1010 to 8.31 × 109 M-1 s-1 at the temperature from 250 to 400 K, respectively, while the ones in the aqueous phase are close to diffusion-controlled rates, with diffusion-corrected rate constants being 1.31 × 109 to 1.27 × 109 M-1 s-1. The formal hydrogen transfer mechanism is the most dominant in the gas phase, whereas the HO•-addition is the most favorable in the aqueous phase. The H-abstraction at two methyl groups and the HO•-addition to C11 and C12 atoms (pyrazole ring), C16 and C18 atoms (benzyl ring) are significant. The short lifetime in the environment, equal to only 4.16 h, requires more attention to this herbicide compound, whereas its lifetime in the aqueous condition varies sharply from half second to several thousand days depending on the HO• concentration. The ecotoxicity estimation of MTZ and its principal transformation products to aquatic organisms suggests that they are harmful or toxic substances. Moreover, the MTZ is a developmental toxicant and mutagenicity-positive, while its decomposed products are developmental toxicants with no mutagenic toxicity. Their bioaccumulation in aquatic organisms is negligible.

New insights into the competition between antioxidant activities and pro-oxidant risks of rosmarinic acid.

Direct and indirect antioxidant activities of rosmarinic acid (RA) based on HOO˙/CH3OO˙ radical scavenging and Fe(iii)/Fe(ii) ion chelation were theoretically studied using density functional theory at the M05-2X/6-311++G(2df,2p) level of theory. First, four antioxidant mechanisms including hydrogen atom transfer (HAT), radical adduct formation (RAF), proton loss (PL) and single electron transfer (SET) were investigated in water and pentyl ethanoate (PEA) phases. Regarding the free radical scavenging mechanism, HAT plays a decisive role with overall rate coefficients of 1.84 × 103 M-1 s-1 (HOO˙) and 4.49 × 103 M-1 s-1 (CH3OO˙) in water. In contrast to PL, RAF and especially SET processes, the HAT reaction in PEA is slightly more favorable than that in water. Second, the [Fe(iii)(H2O)6]3+ and [Fe(ii)(H2O)6]2+ ion chelating processes in an aqueous phase are both favorable and spontaneous especially at the O5, site-1, and site-2 positions with large negative Δr G 0 values and great formation constant K f. Finally, the pro-oxidant risk of RA- was also considered via the Fe(iii)-to-Fe(ii) complex reduction process, which may initiate Fenton-like reactions forming reactive HO˙ radicals. As a result, RA- does not enhance the reduction process when ascorbate anions are present as reducing agents, whereas the pro-oxidant risk becomes remarkable when superoxide anions are found. The results encourage further attempts to verify the speculation using more powerful research implementations of the antioxidant activities of rosmarinic acid in relationship with its possible pro-oxidant risks.

Physicochemical assessment of anionic dye adsorption on bone char using a multilayer statistical physics model.

The statistical physics modeling is a reliable approach to interpret and understand the adsorption mechanism of both organic and inorganic adsorbates. Herein, a theoretical study of the adsorption mechanism of anionic dyes, namely reactive blue 4 (RB4), acid blue 74 (AB74), and acid blue 25 (AB25), on bone char was performed with a multilayer statistical physics model. This model was applied to fit the equilibrium adsorption data of these dyes at 298-313 K and pH 4. Results indicated that the global number of formed dye layers on the bone char varied from 1.62 to 2.24 for RB4, AB74, and AB25 dyes depending on the solution temperature where the saturation adsorption capacities ranged from 0.08 to 0.12 mmol/g. Dye molecular aggregation was also identified for these dyes where dimers and trimers prevailed at different operating conditions especially for adsorbates RB4 and AB74. Adsorption mechanism of these dyes was multimolecular and endothermic with adsorption energies from 10.6 to 20.8 kJ/mol where van der Waals interactions and hydrogen bonding could be present. This investigation contributes to understand the physicochemical variables associated to dye adsorption using low-cost adsorbents as bone char.

Antioxidant and UV-radiation absorption activity of aaptamine derivatives - potential application for natural organic sunscreens.

Antioxidant and UV absorption activities of three aaptamine derivatives including piperidine[3,2-b]demethyl(oxy)aaptamine (C1), 9-amino-2-ethoxy-8-methoxy-3H-benzo[de][1,6]naphthyridine-3-one (C2), and 2-(sec-butyl)-7,8-dimethoxybenzo[de]imidazo[4,5,1-ij][1,6]-naphthyridin-10(9H)-one (C3) were theoretically studied by density functional theory (DFT). Direct antioxidant activities of C1-C3 were firstly evaluated via their intrinsic thermochemical properties and the radical scavenging activity of the potential antioxidants with the HOO˙/HO˙ radicals via four mechanisms, including: hydrogen atom transfer (HAT), single electron transfer (SET), proton loss (PL) and radical adduct formation (RAF). Kinetic calculation reveals that HOO˙ scavenging in water occurs via HAT mechanism with C1 (k app, 7.13 × 106 M-1 s-1) while RAF is more dominant with C2 (k app, 1.40 × 105 M-1 s-1) and C3 (k app, 2.90 × 105 M-1 s-1). Antioxidant activity of aaptamine derivatives can be classified as C1 > C3 > C2. Indirect antioxidant properties based on Cu(i) and Cu(ii) ions chelating activity were also investigated in aqueous phase. All three studied compounds show spontaneous and favorable Cu(i) ion chelating activity with ΔG 0 being -15.4, -13.7, and -15.7 kcal mol-1, whereas ΔG 0 for Cu(ii) chelation are -10.4, -10.8, and -2.2 kcal mol-1 for C1, C2 and C3, respectively. In addition, all compounds show UVA and UVB absorption; in which the excitations are determined mostly as π-π* transition. Overall, the results suggest the potential applications of the aaptamines in pharmaceutics and cosmetics, i.e. as a sunscreen and antioxidant ingredient.

Theoretical Investigation of the Reaction of Pyrene Formation from Fluoranthene.

This paper has investigated the reaction process concerning pyrene formation from fluoranthene in their electronic ground states. Both aromatic compounds are considered as direct soot precursors. The geometrical parameters, the vibrational frequencies, and the zero-point energies have been calculated using the BMK (Boese-Martin for kinetics) method and the 6-311++G(d,p) basis set. More accurate single-point energies have been obtained using BMK/6-311++G(3df,2p) to retrieve thermodynamic properties (ΔrH°(T) and ΔrG°(T)) over a wide temperature range (298-2500 K). The isomerization reaction of fluoranthene to pyrene is exothermic and spontaneous in standard conditions. The transition states and the possible intermediate species have been located on the singlet potential energy surface in order to determine the reaction mechanism. Two different reaction channels have been investigated and characterized by entrance reaction barriers of about 419 and 771 kJ mol-1 for the first and the second reaction pathways, respectively. The present work demonstrates that the first reaction channel is the most energetically favored pathway at high temperatures. Therefore, the kinetic parameters of the forward and reverse first step reactions have been determined in sooting flame conditions.

A CASPT2 theoretical study of the kinetics of the 2-, 3-, and 4-methylbenzylperoxy radical isomerization.

The rate constants of the 2-, 3-, and 4-methylbenzylperoxy isomerization reactions have been computed using the elaborated CASPT2 method. Geometry optimizations and vibrational frequency calculations are performed with two methods (B3LYP and MPW1K) combined with the cc-pVDZ and 6-31+G(d,p) basis sets, respectively. Single-point energy calculations are performed at the CASPT2/ANO-L-VDZP//B3LYP/cc-pVDZ level of theory as recommended by Canneaux et al. (J. Phys. Chem. A 2008, 112, 6045). Canonical transition-state theory with a simple Wigner tunneling correction is used to predict the high-pressure limit rate constants as a function of temperature. They are given by the following relations for the 2-, 3-, and 4-methylbenzylperoxy (MBP) (1,3s) isomerizations and for the 2-methylbenzylperoxy (1,6p) isomerization, respectively: k(2-MBP(1,3s))(600-2000 K) (in s(-1)) = (3.33 x 10(10))T(0.79) exp((-142.6 in kJ mol(-1))/RT); k(3-MBP(1,3s))(600-2000 K) (in s(-1)) = (0.74 x 10(10))T(0.79) exp((-130.7 in kJ mol(-1))/RT); k(4-MBP(1,3s))(600-2000 K) (in s(-1)) = (1.12 x 10(10))T(0.79) exp((-133.6 in kJ mol(-1))/RT); k(2-MBP(1,6p))(600-2000 K) (in s(-1)) = (5.10 x 10(8))T(0.85) exp((-87.1 in kJ mol(-1))/RT). These parameters can be used in the thermokinetic models involving aromatic compounds at high pressure. In the case of the 2-methylbenzylperoxy radical, the (1,6p) H-atom transfer reaction is consistently the most important channel over the studied temperature range, and the (1,3s) H-atom transfer reaction is not energetically favored.

A theoretical study of the kinetics of the benzylperoxy radical isomerization.

The rate constant of the benzylperoxy isomerization reaction has been computed using 54 different levels of theory and has been compared to the experimental value reported at 773 K. The aim of this methodology work is to demonstrate that standard theoretical methods are not adequate to obtain quantitative rate constants for the reaction under study. The use of the elaborated CASPT2 method is essential to estimate a quantitative rate constant. Geometry optimizations and vibrational frequency calculations are performed using three different methods (B3LYP, MPW1K, and MP2) and six different basis sets (6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d,p), 6-311+G(d,p), and cc-pVDZ). Single-point energy calculations are performed with the highly correlated ab initio coupled cluster method in the space of single, double, and triple (pertubatively) electron excitations CCSD(T) using the 6-31G(d,p) basis set, and with the CASPT2 level of theory with the ANO-L-VDZP basis set. Canonical transition-state theory with a simple Wigner tunneling correction is used to predict the high-pressure limit rate constants as a function of temperature. We recommend the use of the CASPT2/ANO-L-VDZP//B3LYP/cc-pVDZ level of theory to compute the temperature dependence of the rate constant of the four-center isomerization of the benzylperoxy radical. It is given by the following relation: k(600-2000 K) (in s (-1)) = (1.29 x 10 (10)) T (0.79) exp[(-133.1 in kJ mol (-1))/ RT]. These parameters can be used in the thermokinetic models involving aromatic compounds at high pressure. This computational procedure can be extended to predict rate constants for other similar reactions where no available experimental data exist.